start

xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0068
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0068
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0067
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0066
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0065
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0063
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0060
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0057
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0053
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0050
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0047
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0046
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0046
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0046
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0047
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0048
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0050
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0051
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0051
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0051
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0012
   -0.0892
    0.9857


xR =

    0.0050
   -0.5079
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0012
   -0.0892
    0.9857


xR =

    0.0049
   -0.5079
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0012
   -0.0892
    0.9857


xR =

    0.0049
   -0.5079
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0012
   -0.0892
    0.9857


xR =

    0.0050
   -0.5079
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0053
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0057
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0061
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0065
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0069
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9856


xR =

    0.0072
   -0.5078
   -0.0081


xk =

    0.0004
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0012
   -0.0892
    0.9856


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0004
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0012
   -0.0892
    0.9856


xR =

    0.0069
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9856


xR =

    0.0069
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9856


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9856


xR =

    0.0073
   -0.5078
   -0.0081


xk =

    0.0004
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9856


xR =

    0.0075
   -0.5078
   -0.0081


xk =

    0.0004
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9856


xR =

    0.0076
   -0.5078
   -0.0081


xk =

    0.0004
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9856


xR =

    0.0077
   -0.5078
   -0.0081


xk =

    0.0004
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9856


xR =

    0.0077
   -0.5078
   -0.0081


xk =

    0.0004
    0.9743
   -0.0001

mean growth 1316

mean_growth_1316 =

    0.1266

mean growth 0712

mean_growth_0712 =

    0.0050

change in potential growth from 2012Q4 to 2016Q3

ans =

    0.0354

average GDP growth 2013 -- 2016 minus 2007 -- 2012 

ans =

    0.1216

GDP growth 2008

Yg2008 =

    0.0111

GDP growth 2009

Yg2009 =

   -0.1085

GDP growth 2010

Yg2010 =

   -0.0522

average GDP growth 2008 -- 2010

ans =

   -0.0499

severity of recession

ans =

   -0.1029

average GDP growth 2010 -- 2012

ans =

    0.0248

(mean(y_tg(1:12)) - log(gamma))*4

ans =

   -0.0154

average revision: (mean(y_tg(12)) - log(gamma))*4

ans =

   -0.0368

GY2009

GY2009 =

    0.1883

GY2010

GY2010 =

    0.1870

GY2011

GY2011 =

    0.1874

GY2012

GY2012 =

    0.1751

change in potential growth from 2009Q4 to 2010Q4

ans =

   -0.0170

change in G

ans =

   -0.1023

GY2010-GY2009

ans =

   -0.0013

GY2011-GY2009

ans =

  -8.8988e-04

GY2012-GY2009

ans =

   -0.0132

RE model@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

xw_RE =

    0.0094
    0.2857
    0.0046

ini_con_RE

ini_con_RE =

   -0.0042    0.0041    0.0038         0

ini_con_RE

ini_con_RE =

   -0.0042    0.0081    0.0071         0

ini_con_RE

ini_con_RE =

   -0.0042    0.0120    0.0101         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0157    0.0126         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0270    0.0218         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0380    0.0298         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0487    0.0368         0

ini_con_RE

ini_con_RE =

         0    0.0592    0.0426         0

ini_con_RE

ini_con_RE =

         0    0.0576    0.0368         0

ini_con_RE

ini_con_RE =

         0    0.0562    0.0311         0

ini_con_RE

ini_con_RE =

         0    0.0547    0.0256         0

ini_con_RE

ini_con_RE =

         0    0.0533    0.0204         0

ini_con_RE

ini_con_RE =

         0    0.0519    0.0153         0

ini_con_RE

ini_con_RE =

         0    0.0506    0.0105         0

ini_con_RE

ini_con_RE =

         0    0.0493    0.0058         0

ini_con_RE

ini_con_RE =

         0    0.0480    0.0013         0

ini_con_RE

ini_con_RE =

         0    0.0468   -0.0030         0

ini_con_RE

ini_con_RE =

         0    0.0456   -0.0071         0

ini_con_RE

ini_con_RE =

         0    0.0444   -0.0111         0

ini_con_RE

ini_con_RE =

         0    0.0433   -0.0149         0

ini_con_RE

ini_con_RE =

         0    0.0422   -0.0185         0

ini_con_RE

ini_con_RE =

         0    0.0411   -0.0220         0

ini_con_RE

ini_con_RE =

         0    0.0401   -0.0254         0

ini_con_RE

ini_con_RE =

         0    0.0390   -0.0286         0

ini_con_RE

ini_con_RE =

         0    0.0380   -0.0316         0

ini_con_RE

ini_con_RE =

         0    0.0371   -0.0346         0

ini_con_RE

ini_con_RE =

         0    0.0361   -0.0374         0

ini_con_RE

ini_con_RE =

         0    0.0352   -0.0401         0

ini_con_RE

ini_con_RE =

         0    0.0343   -0.0426         0

ini_con_RE

ini_con_RE =

         0    0.0334   -0.0451         0

ini_con_RE

ini_con_RE =

         0    0.0326   -0.0474         0

ini_con_RE

ini_con_RE =

         0    0.0317   -0.0496         0

ini_con_RE

ini_con_RE =

         0    0.0309   -0.0517         0

ini_con_RE

ini_con_RE =

         0    0.0302   -0.0538         0

ini_con_RE

ini_con_RE =

         0    0.0294   -0.0557         0

ini_con_RE

ini_con_RE =

         0    0.0286   -0.0575         0

ini_con_RE

ini_con_RE =

         0    0.0279   -0.0592         0

ini_con_RE

ini_con_RE =

         0    0.0272   -0.0609         0

ini_con_RE

ini_con_RE =

         0    0.0265   -0.0624         0


hp2 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [1 0 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'o'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×56 double]
              YData: [1×56 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('hp2', 'var'), matlab.graphics.internal.getForDisplay('hp2', hp2, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('hp2'), end">all properties</a>


ythand = 

  <a href="matlab:helpPopup matlab.graphics.primitive.Text" style="font-weight:bold">Text</a> with properties:

                 String: ''
               FontSize: 8.8000
             FontWeight: 'normal'
               FontName: 'Helvetica'
                  Color: [0.1500 0.1500 0.1500]
    HorizontalAlignment: 'center'
               Position: [-15.5341 7.5000 -1]
                  Units: 'data'

  Show <a href="matlab:if exist('ythand', 'var'), matlab.graphics.internal.getForDisplay('ythand', ythand, 'matlab.graphics.primitive.Text'), else, matlab.graphics.internal.getForDisplay('ythand'), end">all properties</a>


ythand = 

  <a href="matlab:helpPopup matlab.graphics.primitive.Text" style="font-weight:bold">Text</a> with properties:

                 String: ''
               FontSize: 8.8000
             FontWeight: 'normal'
               FontName: 'Helvetica'
                  Color: [0.1500 0.1500 0.1500]
    HorizontalAlignment: 'center'
               Position: [-15.4512 -3.7547 -1.0000]
                  Units: 'data'

  Show <a href="matlab:if exist('ythand', 'var'), matlab.graphics.internal.getForDisplay('ythand', ythand, 'matlab.graphics.primitive.Text'), else, matlab.graphics.internal.getForDisplay('ythand'), end">all properties</a>

change in output gaps

ans =

   -0.0244

change in debt GDP ratio

ans =

    0.0342

grow_err(6)

ans =

   -0.0244

CAB 2009Q4

ans =

    0.0060

change in CAB: 2009Q4 -- 2012Q4

ans =

    0.0036

change in potential output growth: 2009Q4 -- 2012Q4

ans =

   -0.4832

growth_fore_err

growth_fore_err =

   -0.0762


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0068
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0068
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0067
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0066
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0065
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0063
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0060
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0057
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0053
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0050
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0048
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0047
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0047
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0047
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0047
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0095
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0049
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0050
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0052
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0053
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0055
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0057
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0094
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0059
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0061
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0063
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0064
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0066
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0067
   -0.5078
   -0.0082


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0068
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0093
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0069
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0070
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0070
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0013
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001


xw =

    0.0092
    0.2857
    0.0046


xb =

    0.0014
   -0.0892
    0.9857


xR =

    0.0071
   -0.5078
   -0.0081


xk =

    0.0003
    0.9743
   -0.0001

mean growth 1316

mean_growth_1316 =

    0.0958

mean growth 0712

mean_growth_0712 =

    0.0203

change in potential growth from 2012Q4 to 2016Q3

ans =

    0.0266

average GDP growth 2013 -- 2016 minus 2007 -- 2012 

ans =

    0.0754

GDP growth 2008

Yg2008 =

    0.0111

GDP growth 2009

Yg2009 =

   -0.1035

GDP growth 2010

Yg2010 =

   -0.0503

average GDP growth 2008 -- 2010

ans =

   -0.0476

severity of recession

ans =

   -0.1006

average GDP growth 2010 -- 2012

ans =

    0.0538

(mean(y_tg(1:12)) - log(gamma))*4

ans =

   -0.0150

average revision: (mean(y_tg(12)) - log(gamma))*4

ans =

   -0.0357

GY2009

GY2009 =

    0.1890

GY2010

GY2010 =

    0.1877

GY2011

GY2011 =

    0.1874

GY2012

GY2012 =

    0.1876

change in potential growth from 2009Q4 to 2010Q4

ans =

   -0.0161

change in G

ans =

   -0.0956

GY2010-GY2009

ans =

   -0.0012

GY2011-GY2009

ans =

   -0.0016

GY2012-GY2009

ans =

   -0.0014

RE model@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

xw_RE =

    0.0094
    0.2857
    0.0046

ini_con_RE

ini_con_RE =

   -0.0042    0.0041    0.0038         0

ini_con_RE

ini_con_RE =

   -0.0042    0.0081    0.0071         0

ini_con_RE

ini_con_RE =

   -0.0042    0.0120    0.0101         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0157    0.0126         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0270    0.0218         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0380    0.0298         0

ini_con_RE

ini_con_RE =

   -0.0120    0.0487    0.0368         0

ini_con_RE

ini_con_RE =

         0    0.0592    0.0426         0

ini_con_RE

ini_con_RE =

         0    0.0576    0.0368         0

ini_con_RE

ini_con_RE =

         0    0.0562    0.0311         0

ini_con_RE

ini_con_RE =

         0    0.0547    0.0256         0

ini_con_RE

ini_con_RE =

         0    0.0533    0.0204         0

ini_con_RE

ini_con_RE =

         0    0.0519    0.0153         0

ini_con_RE

ini_con_RE =

         0    0.0506    0.0105         0

ini_con_RE

ini_con_RE =

         0    0.0493    0.0058         0

ini_con_RE

ini_con_RE =

         0    0.0480    0.0013         0

ini_con_RE

ini_con_RE =

         0    0.0468   -0.0030         0

ini_con_RE

ini_con_RE =

         0    0.0456   -0.0071         0

ini_con_RE

ini_con_RE =

         0    0.0444   -0.0111         0

ini_con_RE

ini_con_RE =

         0    0.0433   -0.0149         0

ini_con_RE

ini_con_RE =

         0    0.0422   -0.0185         0

ini_con_RE

ini_con_RE =

         0    0.0411   -0.0220         0

ini_con_RE

ini_con_RE =

         0    0.0401   -0.0254         0

ini_con_RE

ini_con_RE =

         0    0.0390   -0.0286         0

ini_con_RE

ini_con_RE =

         0    0.0380   -0.0316         0

ini_con_RE

ini_con_RE =

         0    0.0371   -0.0346         0

ini_con_RE

ini_con_RE =

         0    0.0361   -0.0374         0

ini_con_RE

ini_con_RE =

         0    0.0352   -0.0401         0

ini_con_RE

ini_con_RE =

         0    0.0343   -0.0426         0

ini_con_RE

ini_con_RE =

         0    0.0334   -0.0451         0

ini_con_RE

ini_con_RE =

         0    0.0326   -0.0474         0

ini_con_RE

ini_con_RE =

         0    0.0317   -0.0496         0

ini_con_RE

ini_con_RE =

         0    0.0309   -0.0517         0

ini_con_RE

ini_con_RE =

         0    0.0302   -0.0538         0

ini_con_RE

ini_con_RE =

         0    0.0294   -0.0557         0

ini_con_RE

ini_con_RE =

         0    0.0286   -0.0575         0

ini_con_RE

ini_con_RE =

         0    0.0279   -0.0592         0

ini_con_RE

ini_con_RE =

         0    0.0272   -0.0609         0

ini_con_RE

ini_con_RE =

         0    0.0265   -0.0624         0

+++++++++++++++++++multiplier++++++++++++++++

multiplier =

  Columns 1 through 13

         0         0         0         0         0         0         0         0    2.3866    2.9106    3.4096    3.6780    3.8375

  Columns 14 through 26

    3.9357    3.9946    3.9560    2.9009    2.3287    2.2407    2.3257    2.4544    2.6052    2.7692    2.9402    3.1156    3.2388

  Columns 27 through 39

    3.3128    3.3506    3.3432    3.2346    3.0911    2.9861    2.9542    2.9962    3.0087    2.9963    2.9619    2.9072    2.8333

  Column 40

    2.7408


multiplier1 =

  Columns 1 through 13

         0         0         0         0         0         0         0         0    4.7286    4.8138    4.8799    4.9111    4.9278

  Columns 14 through 26

    4.9369    4.9415    4.9402    4.8757    4.7642    4.6333    4.5041    4.3876    4.3006    4.2518    4.2412    4.2639    4.2834

  Columns 27 through 39

    4.2931    4.2948    4.2809    4.2176    4.1262    4.0411    3.9903    3.9883    3.9794    3.9645    3.9444    3.9199    3.8917

  Column 40

    3.8606


multiplier_RE =

  Columns 1 through 21

     0     0     0     0     0     0     0     0   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN

  Columns 22 through 40

   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN   NaN


multiplier_RE1 =

  Columns 1 through 13

         0         0         0         0         0         0         0         0   -3.5958   -4.3571   -5.3940   -6.8889   -9.2309

  Columns 14 through 26

  -13.4243  -23.0980  -69.2396   81.2048   26.8850   16.5594   12.1864    9.7705    8.2383    7.1799    6.4050    5.8133    5.3467

  Columns 27 through 39

    4.9693    4.6578    4.3964    4.1738    3.9822    3.8153    3.6688    3.5391    3.4236    3.3199    3.2265    3.1418    3.0647

  Column 40

    2.9942


years =

        2008        2009        2010        2011        2012        2013        2014        2015        2016        2017


h1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 1 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [733682 734047 734412 734777 735143 735508]
              YData: [-7.4624 -18.4008 -23.3022 -26.0141 -31.1630 -35.0287]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h1', 'var'), matlab.graphics.internal.getForDisplay('h1', h1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h1'), end">all properties</a>


h1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h1', 'var'), matlab.graphics.internal.getForDisplay('h1', h1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h1'), end">all properties</a>


y_with =

  Columns 1 through 13

   -0.4190   -1.3173   -2.6033   -4.1904   -6.7764  -10.4170  -15.0112  -20.3433  -24.9215  -28.1331  -30.0043  -30.8607  -30.9298

  Columns 14 through 26

  -30.3377  -29.2073  -27.6891  -26.6927  -26.7254  -27.5643  -28.7966  -30.0438  -30.6606  -29.8571  -27.1068  -22.3057  -16.8207

  Columns 27 through 39

  -11.7781   -7.3942   -3.9863   -2.9592   -4.2458   -6.1270   -6.6898   -4.5576   -1.6652    0.6588    2.4467    3.7395    4.5845

  Column 40

    5.0329


ylhand = 

  <a href="matlab:helpPopup matlab.graphics.primitive.Text" style="font-weight:bold">Text</a> with properties:

                 String: ''
               FontSize: 8.8000
             FontWeight: 'normal'
               FontName: 'Helvetica'
                  Color: [0.1500 0.1500 0.1500]
    HorizontalAlignment: 'center'
               Position: [7.3318e+05 -20.0000 -1.0000]
                  Units: 'data'

  Show <a href="matlab:if exist('ylhand', 'var'), matlab.graphics.internal.getForDisplay('ylhand', ylhand, 'matlab.graphics.primitive.Text'), else, matlab.graphics.internal.getForDisplay('ylhand'), end">all properties</a>

average output loss over three years

ans =

   -1.9208

cumulative output loss over six years

ans =

  -28.5756


h2 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h2', 'var'), matlab.graphics.internal.getForDisplay('h2', h2, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h2'), end">all properties</a>


h3 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '--'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h3', 'var'), matlab.graphics.internal.getForDisplay('h3', h3, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h3'), end">all properties</a>


h3 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h3', 'var'), matlab.graphics.internal.getForDisplay('h3', h3, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h3'), end">all properties</a>


h4 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '--'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h4', 'var'), matlab.graphics.internal.getForDisplay('h4', h4, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h4'), end">all properties</a>


h4 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h4', 'var'), matlab.graphics.internal.getForDisplay('h4', h4, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h4'), end">all properties</a>


h5 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '--'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h5', 'var'), matlab.graphics.internal.getForDisplay('h5', h5, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h5'), end">all properties</a>


hi1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('hi1', 'var'), matlab.graphics.internal.getForDisplay('hi1', hi1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('hi1'), end">all properties</a>


hi1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '--'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('hi1', 'var'), matlab.graphics.internal.getForDisplay('hi1', hi1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('hi1'), end">all properties</a>


h1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 1 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [734412 734777 735143 735508]
              YData: [0 0 0 0]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h1', 'var'), matlab.graphics.internal.getForDisplay('h1', h1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h1'), end">all properties</a>

g_with1

g_with1 =

  Columns 1 through 13

   -0.0036   -0.0108   -0.0211   -0.0338   -0.0549   -0.0949   -0.2221   -0.3255   -0.3961   -0.3151   -0.2381   -0.2494   -0.2554

  Columns 14 through 26

   -0.2571   -0.2554   -0.2725   -0.7877   -1.3645   -1.8541   -2.2115   -2.4498   -2.3552   -1.8880   -1.1428   -0.2572   -0.0520

  Columns 27 through 39

   -0.0261   -0.0048   -0.1844   -1.1869   -1.8904   -1.9188   -1.2074    0.0693    0.0822    0.0914    0.0972    0.1000    0.1001

  Column 40

    0.0980


h2 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h2', 'var'), matlab.graphics.internal.getForDisplay('h2', h2, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h2'), end">all properties</a>


h2 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '--'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h2', 'var'), matlab.graphics.internal.getForDisplay('h2', h2, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h2'), end">all properties</a>


h1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 1 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [734047 734412 734777 735143 735508 735873 736238 736604]
              YData: [0 0 0 0 0 0 0 0]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h1', 'var'), matlab.graphics.internal.getForDisplay('h1', h1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h1'), end">all properties</a>


hp = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 0 1]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'o'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('hp', 'var'), matlab.graphics.internal.getForDisplay('hp', hp, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('hp'), end">all properties</a>


hpp = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [1 0 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'diamond'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('hpp', 'var'), matlab.graphics.internal.getForDisplay('hpp', hpp, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('hpp'), end">all properties</a>


y2hand = 

  <a href="matlab:helpPopup matlab.graphics.primitive.Text" style="font-weight:bold">Text</a> with properties:

                 String: ''
               FontSize: 9.9000
             FontWeight: 'normal'
               FontName: 'Helvetica'
                  Color: [0.1500 0.1500 0.1500]
    HorizontalAlignment: 'center'
               Position: [7.3294e+05 3.8147e-05 -1]
                  Units: 'data'

  Show <a href="matlab:if exist('y2hand', 'var'), matlab.graphics.internal.getForDisplay('y2hand', y2hand, 'matlab.graphics.primitive.Text'), else, matlab.graphics.internal.getForDisplay('y2hand'), end">all properties</a>


hpp1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [1 0 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'diamond'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×40 double]
              YData: [1×40 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('hpp1', 'var'), matlab.graphics.internal.getForDisplay('hpp1', hpp1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('hpp1'), end">all properties</a>


y3hand = 

  <a href="matlab:helpPopup matlab.graphics.primitive.Text" style="font-weight:bold">Text</a> with properties:

                 String: ''
               FontSize: 9.9000
             FontWeight: 'normal'
               FontName: 'Helvetica'
                  Color: [0.1500 0.1500 0.1500]
    HorizontalAlignment: 'center'
               Position: [7.3294e+05 1.9073e-05 -1.0000]
                  Units: 'data'

  Show <a href="matlab:if exist('y3hand', 'var'), matlab.graphics.internal.getForDisplay('y3hand', y3hand, 'matlab.graphics.primitive.Text'), else, matlab.graphics.internal.getForDisplay('y3hand'), end">all properties</a>


h1 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [0 1 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'none'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [733682 734047 734412 734777 735143 735508]
              YData: [0 0 0 0 0 0]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('h1', 'var'), matlab.graphics.internal.getForDisplay('h1', h1, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('h1'), end">all properties</a>


hp2 = 

  <a href="matlab:helpPopup matlab.graphics.chart.primitive.Line" style="font-weight:bold">Line</a> with properties:

              Color: [1 0 0]
          LineStyle: '-'
          LineWidth: 0.5000
             Marker: 'o'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×56 double]
              YData: [1×56 double]
              ZData: [1×0 double]

  Show <a href="matlab:if exist('hp2', 'var'), matlab.graphics.internal.getForDisplay('hp2', hp2, 'matlab.graphics.chart.primitive.Line'), else, matlab.graphics.internal.getForDisplay('hp2'), end">all properties</a>


ythand = 

  <a href="matlab:helpPopup matlab.graphics.primitive.Text" style="font-weight:bold">Text</a> with properties:

                 String: ''
               FontSize: 8.8000
             FontWeight: 'normal'
               FontName: 'Helvetica'
                  Color: [0.1500 0.1500 0.1500]
    HorizontalAlignment: 'center'
               Position: [-15.5341 7.5000 -1]
                  Units: 'data'

  Show <a href="matlab:if exist('ythand', 'var'), matlab.graphics.internal.getForDisplay('ythand', ythand, 'matlab.graphics.primitive.Text'), else, matlab.graphics.internal.getForDisplay('ythand'), end">all properties</a>


ythand = 

  <a href="matlab:helpPopup matlab.graphics.primitive.Text" style="font-weight:bold">Text</a> with properties:

                 String: ''
               FontSize: 8.8000
             FontWeight: 'normal'
               FontName: 'Helvetica'
                  Color: [0.1500 0.1500 0.1500]
    HorizontalAlignment: 'center'
               Position: [-15.4512 -3.6641 -1]
                  Units: 'data'

  Show <a href="matlab:if exist('ythand', 'var'), matlab.graphics.internal.getForDisplay('ythand', ythand, 'matlab.graphics.primitive.Text'), else, matlab.graphics.internal.getForDisplay('ythand'), end">all properties</a>

change in output gaps

ans =

   -0.0242

change in debt GDP ratio

ans =

    0.0343

grow_err(6)

ans =

   -0.0244

CAB 2009Q4

ans =

    0.0060

change in CAB: 2009Q4 -- 2012Q4

ans =

    0.0036

change in potential output growth: 2009Q4 -- 2012Q4

ans =

   -0.4832

growth_fore_err

growth_fore_err =

    0.0107

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
RE solution
ordering: rk wage output mp hours cons invst caput cap gama
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
